Comment on "High— order contamination in the tail of gravitational collapse 
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We confront the predictions of S. Hod, Phys. Rev. D 60, 104053 (1999) for tiie late-time decay 
rate of black hole perturbations with numerical data. Specifically, we ask two questions: First, are 
corrections to the Price tail dominated by logarithmic terms, as predicted by Hod? Second, if there 
were logarithmic correction terms, do they take the specific form predicted in Hod's paper? The 
answer to both questions is "no." 
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The late-time tail of black hole perturbations has been 
studied in detail HH Hills,, 6, 7,, 8]. The leading-order 
(in inverse-time) behavior of the tail is well understood: 
it is given by a power-law in inverse-time, the power in- 
dex depending on the multipole order of the initial per- 
turbation. Specifically, the late-time behavior is given 
asymptotically by i^™. For initial data with compact 
support, m — 2t + i, i being the multipole order. (For 
a rotating black hole mode coupling can generate lower 
values of m if not disallowed 6].) In Ref. Ji], the ques- 
tion of higher-order corrections is addressed. Specifically, 
Ref. IqI argues that the correction is dominated by a log- 
arithmic term. Indeed, naive expectations support this 
prediction: as the effective potential that scatters the 
waves includes a logarithmic term (when expressed in 
"natural" coordinates, such as the Regge-Wheeler "tor- 
toise" coordinate) , it is natural to expect that logarithmic 
terms would appear also in the scattered waves. Notwith- 
standing these expectations, in this Comment we show 
that the correction terms show no evidence of a logarith- 
mic dependence, and that it is dominated by a simple 
power series. In fact, Ching et al Q showed that the 
Schwarzschild potential is exceptional, in that the lead- 
ing behavior is not logarithmic, although the potential 
is. Most logarithmic potentials would induce logarith- 
mic tails. In this Comment we argue, and demonstrate 
numerically, that the conclusions of Ching et al apply 
not just to the leading order tail, but also to the leading 
subdominant term. 

That simple higher-order power-series terms exist in 
the tail is a trivial observation: It is well known that 
the leading order term in an inverse-time expansion is 
given by A small shift in the origin of time, i.e., the 
transformation t ^ t — t^^ {t ^ t^^, as we are interested in 
the late-time limit, as t 00) will result in an asymp- 
totic behavior given by {t ~ to)-'"' - i"™ + mto ^"("+1) . 
Namely, a simple shift in the origin of time gives rise to a 
simple power series expansion of the late-time tail, that 
has nothing to do with black hole scattering. In actual 



simulations, the origin of time is typically chosen arbi- 
trarily. That is, as no special attention is typically given 
to the question of where to place the origin of time, one 
would expect generic simulations to always include such 
simple-series terms. The interesting question then is not 
whether simple-series terms exist, but rather whether 
they are dominated by logarithmic terms, as predicted 
by Ref. 9J. We expect that the answer to the latter 
question is "no." These expectations are supported by a 
failure to reproduce logarithmic terms for the late-time 
tail, both by the present authors and by others p^. 

In the picture that we present in this Comment, the 
effect of black-hole scattering would then be to add to 
the correction term, i.e., the behavior of the late-time 
field would now be given by i"™ + {a + mto)t~^"''^^\ 
where a is a constant (that may depend on the location 
of the evaluation point and the parameters of the pulse). 
Distinguishing between a pure power series and a power 
series that includes logarithmic terms is not easy. Specif- 
ically, it is not easy to find from the behavior of the field 
itself whether logarithmic correction terms appear or not, 
for the simple reason that the function In t increases only 
very slowly with t, and there are practical limits to how 
long one can evolve the data numerically. Instead, it is 
useful to calculate auxiliary quantities — some of which 
have been introduced before — and look for such quanti- 
ties for which the signature of the logarithmic correction 
terms is easily recognizable. In this Comment we fo- 
cus on three such quantities: (a) the local power index 
n := -tiP/^ "S], (b) S := (n - Uoo) (t/M)/ lii{t/M) 9^, 
and (c) A :— {n — n,^) (t/M), where :— n{t 00). 
Here we assume, without loss of generality, that the per- 
turbation is that of a Schwarzschild black hole of mass 
M, and that the perturbing field is a spherical {£ — 0) 
scalar field, so as to make rioo = 3. In what follows, we 
assign a subscript "A''" to numerically evaluated data, 
and a subscript "iJ" to quantities predicted by Ref. Q . 

We propose that the late-time behavior is given by a 
local power index n as follows: 
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FIG. 1: The Holder Lp vector norms for the local power index. 
Shown are the Lp norms for p = 1, 2, oo for a vector of ten 
equally spaced data points at r* = between t — QOOOAf 
and t = lOOOOM. Each norm corresponds to the vector of 
differences of the data obtained for a grid density of 2N data 
points per M and those obtained for a grid density of A'^ data 
points per M. Data are shown for A'^ — 1,2,4,8. The Li 
(sum) norm is shown in squares, the L2 (Euclidean) norm in 
circles, and the Loo (maximum) norm in diamonds. 

and the prediction of Ref. |^ is that 

„ , \n(t/M) M 

for some dimensionless constants A, B; Ah, Bh- In fact, 
Ref. Q predicts that Ah = 12 exactly. One way to 
distinguish between the two predictions is to plot n^r 
as a function of M/t, and test whether it approaches a 
straight line asymptotically as M/t — > 0. This is done in 
Fig. 13 which displays the numerical data points and the 
best fit curves based on our ansatz and the ansatz 
of Ref. "^0. The best fit is done using linear regres- 
sion and the least squares method. We did not consider 
here the numerical noise in individual data points. We 
used a 1 + ID code in double-null coordinates whose con- 
vergence is second order. We specified initially incoming 
scalar-field perturbations on the characteristic hypcrsur- 
face, with profile 

^ _ {v - vi){v3 - v) ^ 

[{V2 - Vi){v3 - V2)_ 

for vi < V < V3, and ■0 = otherwise, and with 
V2 = {vi + V3)/2. Here, v is advanced time. In prac- 
tice, we chose vi — lOM and U3 = 30M. Although 
the quantitative results (e.g., the value of the parameters 
A, B, Ah, Bh) depend on the choice of initial data, the 
qualitative behavior — and our conclusions — are indepen- 
dent of the choice of the initial data. The stability and 




FIG. 2: The diagnostic A as a function oi M/t for different 
grid densities (in terms of the number A'^ of grid point per M): 
Bold solid curve (A'^ = 1), dotted curve (A'^ = 2), dash-dotted 
curve (A'^ = 4), dashed curve (A'' = 8), and thin solid curve 
{N = 16). 
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FIG. 3: The local power index n — 3 as a function of M/t. 
The numerical data points are depicted by circles; The solid 
line is a best fit curve based on the ansatz (0, and the two 
dotted curves are the 3 a confidence curves; The dashed curve 
is a best fit curve based on the ansatz |(2J , and the two dash- 
dotted curves are the 3 a confidence curves. 



second-order convergence of the code is demonstrated in 
Fig. n that shows the Holder Lp vector norms of the lo- 
cal power index at various grid densities, and in Fig. [21 
that shows the diagnostic A as a function of time for 
different grid densities. Other indicators behave simi- 
larly. However, all our diagnostics depend crucially on 
the local power index, whose calculation requires the nu- 
merical differentiation of the field. In fact, we found that 
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FIG. 4: (5 as a function of M/t. Left panel: The circles are 
the numerical data points; The dotted curve is a first-order 
best fit curve, and the solid line a second-order best fit curve, 
based on the ansatz JIJ; The dashed curve is a second-order 
best fit curve based on the ansatz 0. Right panel: The 
difference between the three best fit curves and the numerical 
data. 




FIG. 5: ^ — 12 as a function of M/t. The circles are the 
numerical data points, the dashed line is a best fit curve based 
on the ansatz with Ah = 12, and the two dash-dotted 
curves are the 3 a confidence curves. The solid curve is a best 
fit curved based on the ansatz I^Jl with the corresponding 3 a 
confidence curves depicted by the dotted curves. 



the common double precision floating point arithmetic 
results in round off noise that prevents us from finding 
the phenomenon of interest easily. We solve the problem 
by using quadrupole precision floating point arithmetic. 
Our code exhibits clear second-order convergence glob- 
ally, throughout the entire domain of the computation. 



FIG. 6: A as a function of M/t. The circles are the numerical 
data point, the solid curve is a best fit curve based on the 
ansatz @ with the corresponding 3 a confidence curves shown 
as the dotted curves, and the dashed curve is the best fit curve 
based on the ansatz ||7|l, with the corresponding 3 a confidence 
curves shown by the dash-dotted curves. 



Figure |31 displays the local power index as a funtion 
of time. It suggests that our ansatz fits the numeri- 
cal data better than the ansatz of Ref. |^. Indeed, our 
curve has a corresponding squared correlation coefficient 
i?2 = 0.99985, whereas = 0.98589. TableHlshows the 
best fit values for the parameters of the curves. Recall 
that not only does Ref. @ make the claim that there are 
logarithmic correction terms, it also predicts the value 
of the expansion parameter Ah to equal 12. However, 
the best fit finds that Ah = 13.77 ± 0.17, which deviates 
from the prediction of Ref. ^9] by over 10 cr. Although 
our best fit curves agrees with the numerical data to a 
greater extent, it is not easy to distinguish between an 
asymptotically linear curve (in t~^), and a curve that be- 
haves asymptotically like \nt/t. In fact, one can construe 
the definition for n as a differential equation for the field 
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and if one can show that n is asymptotically linear in 
our ansatz is proved, and the prediction of Ref. 
falsified. However, as we have just seen, it is not easy to 
decide on this question from direct observation of n. We 
therefore consider other diagnostics. 

The next diagnostic is 5, which was first introduced in 
Ref. 0. Our prediction is that 

.^^+ ^ (4) 
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TABLE L Determination of the parameters A, B and Ah, Bh, 
and the value of the squared correlation coefhcients and 
Rjj for either ansatz, using the various diagnostics. The data 
above the single horizontal line refers to our ansatz of no 
logarithmic terms^and the data below the line refers to the 
prediction of Ref. (a): In this case, because the predicted 
A is constant, strictly speaking — 0. (b): In this case 
Ah anfi-correlates with the numerical data; the former is 
a decreasing function of M/t, and the latter an increasing 
function. 



Determination 



from 


A or Ah 


B or Bh 


or Rjj 


n — 3 (1st order) 


99.51 ±0.13 




0.999851 


5 (1st order) 


98.99 ±0.12 




0.994955 


5 (2nd order) 


97.76 ± 0.05 


2010 ± 73 


0.999901 


A (1st order) 


98.89 ± 0.13 




(a) 


A (2nd order) 


97.708 ± 0.008 


2059 ± 12 


0.99946 


analytical 


12 






nn -3 (1st order) 13.77 ± 0.17 




0.98589 


5h (1st order) 


1.1 ±0.23 




(a) 


5h (2nd order) 


-0.93 ±0.09 


106.02 ± 0.67 0.999365 


Sh - 12 




8.8 ± 1.7 


0.154808 


Ah (1st order) 


12.96 ± 0.25 




(b) 


Ah (2nd order) 


-0.94 ±0.07 


106.1 ±0.5 


0.928968 



Figure^ shows S (left panel) and 6 — Sn (right panel) as 
functions of M/t, for our first- and second-order predic- 
tions, and the second-order prediction of Ref. 0- As in 
the case of n — 3, our ansatz appears to agree better with 
the numerical data. Indeed, we find from our first-order 
curve that i?2 ^ 0.994955 and from the second-order 
curve = 0.999901. The corresponding values for the 
predictions of Ref. Ig are Rh^O and Rj^ = 0.999365. 
While the first-order prediction of Ref. 9] fails spec- 
tacularly, its second-order prediction fits the numerical 
data quite well. (This, in itself, is not surprising, as the 
second-order term of Ref. [3] is identical with our first- 
order term.) Also, the best fit determination for the pa- 
rameters of the curve does not appear to converge for 
the predictions of Ref. Q : the first-order determination 
is Ah — 1.1 ± 0.2 and the second-order determination 



is Ah = -0.93 ± 0.09. Recall that Ref. ^ predicts not 
only logarithmic terms, but also predicts Ah = 12. The 
latter result for Ah deviates from the value of 12 pre- 
dicted in Ref. ^ by 144 (7. Next, we assume Ah = 12, 
and present in Fig.|5l(5 — 12 as a function of M/t. We find 
Rjf — 0.154808. Clearly, this prediction is not supported 
by the numerical data. 

Lastly, we use A as a diagnostic. Our prediction for A 

is 

A = A + B^ + --- (6) 

whereas 

Ah = Ah Mt/M) + Bh + --- . (7) 

Figure ini displays A as a function of M/t for our second- 
order_prediction, and the second-order prediction of 
Ref. I9l . Notice that the curvature of the curve predicted 
by Ref. ^ has the wrong concavity (cf. also with Fig.[5J). 
We find that R^ = 0.99946, and Rjj = 0.928968. Also, 
Ah = —0.94 ± 0.07, which disagrees with the predicted 
value of 12 by 198.2 ct. 

Our results hold also for the case of a Kerr black hole. 
Ref. has shown that indeed the local power index is 
asymptotically linear in t~^ also for the case of Kerr per- 
turbations. Using Eq. (O, we conclude that the late- 
time field is a simple power series in inverse time also 
for the case of Kerr. (Notice that this question is sep- 
arate from the question of mode couplings in the Kerr 
case.) Our results show that a simple power-series fits 
numerical data much better than an ansatz that include 
logarithmic terms. In addition, determination of the pa- 
rameters using different diagnostics is inconsistent, when 
the prediction of Ref. is assumed. Moreover, the spe- 
cific determination of Ref. that Ah = 12 is shown here 
to be incorrect to a very high confidence level. 
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